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Fig. 1. Geographical location of the study area. 


geographical distribution, ecological and ethnobotanical adaptation of a group of plants with geography, 
ecology, climate and human settlements... ” (Guarino et al. 2002), can be used to analyze this distribution. 
Ecological and geographical information on species is fundamental for conservation, planning, and prediction 
(Ferner 2002; Funk & Richardson 2002; Rushton et al. 2004), as well as to reach an understanding of factors 
which determine the patterns of spatial distribution of biodiversity (Elith et al. 2006; Ricklefs 2004). In gen¬ 
eral, occurrence data for the great majority of species are scarce and of presence-only since they generally come 
rom random samplings that study flora and fauna. Therefore, determining the species distribution using data 
generated in such manner is inadequate for many applications (Cortes & Herrera 2011). Spatial distribution 
stu ies attempt to provide predictions of the distribution of species by connecting their presence with envi¬ 
ronmental factors. These have been applied in the study of relationships between environmental parameters 
species ric ness (Herrera & Cortes 2010; MacNally & Fleishman 2004) and invasion by non-native spe¬ 
cies oosy ). e present work attempts to analyze a large quantity of data compiled through flora and 
launa studies in order to imnmvp _r .... r 


... * -/ viuaiiiliv U 

fauna slud.es m order to improve current knowledge of species distribution. 

To this end Jocatton data (latitude and longitude) for approximately 500 grass species have been obtained 

“r ,y r eS frOIn the S,at “ ° f Chih “ h ^ and Zacatecas (Herrera 2001; Her- 

1 TV 0 " f errera Cortes 2009,2010; Herrera et al. 2010; Cortes & Herrera 2011; Herrera et al 

in nrCD. ). 1 his information urac _ \ __ j % i _ 


orocc enpripe p >u , u lu Ufcrierr nme the climatic requirements tor growtn ui 

duced non-nati™" 1 Th* ° 4 grasses from north central Mexico were chosen: 3 native and 1 intro- 
ecosystems allowinp C fh S & r . dl * tnbutlon data of these 4 species indicate that they occupy divers* 

based on the climatiOa T m J2 ° l * hypothesis that potential distribution of plants can be determined 
based on he c i inane factors of the physical environment. 
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porteri, and Melinis repens ) collected in the study area were analyzed with Geographic Information Systems 
(GIS) software, specifically the DIVA-GIS 7.1.7.2 version (http://www.diva-gis.org/), and GVS1G 1.9 for Vista 
(http://www.gvsig.org/), using geo-statistics maps of the States of Chihuahua, Durango, and Zacatecas (INEG1 
geo-statistics framework) together with Bioclimatic data of 30-second arc resolution from the Worldclim da¬ 
tabase (http://worldclim.org/current). The GIS software packages were developed for the analysis of species 
diversity and distribution in order to clarify geographical, ecological, and genetic standards (Hijmans et al. 
2004) and contain the necessary algorithms to generate the required data in a simple manner. 

The study area encompasses the territory of the States of Chihuahua, Durango, and Zacatecas in the north 
central Mexico, occupying approximately 446,000 km 2 (Fig. 1) located between the extreme coordinates 
21°02'24" to 31°46'48"N and 109°04T2" to 100°44'24"W. 


This area is also located within the physiographical provinces of the Sierra Madre Occidental and Sierras 
del Norte, as well as partially within the Mesa del Centro (INEGI 2003b). As such, the topography of the study 
area is highly varied and rugged with altitude variations between less than 500 m and up to 3200 m. Its latitude 
position, the distance to and influence of the maritime zone, and the topographical characteristics of this zone 
determine strong variation in climatic factors such as temperature and precipitation. Within this region arc 
low-altitude areas with an annual temperature ranging between 24 to 26°C and high-alt itude areas with annual 
temperatures ranging between 8 and 10°C. In terms of precipitation the study area has dry zones with little 
annual precipitation (200 mm) as well as sub-humid zones that have an annual precipitation between 1000 
and 1500 mm. Interrelations between these and other factors cause the study area to have varied plant commu¬ 
nities and soil cover from zones with pine and oak forests to grasslands and xerophilus scrubs (microphyll 
desert scrub, rosette desert scrub, crassicaule scrub), as well as small pockets dominated by halophilous scrub 
communities or halophilous grasslands, due to salty soils (INEGI 2003a). 

The potential distribution analysis was carried out using a subset of four (out of 500) species belonging to 
diverse environments in order to determine the potential distribution and climatic characteristics of the distri¬ 
bution area. A total of three native species were selected, Muhlenbergio pevuvio.no from temperate semiarid 
forests; Boutelouo grocilis from grasslands; Muhlenbergio porteri from the xeric region; as well as one nonnative 
introduction, Melinis repens (Zizka rose Natal grass). The data used to carry out our analysis came from speci¬ 


mens kept at several herbaria in Mexico (CHAPA, CIIDIR, ENCB, HUAA, 1EB, MEXU, SLPM) as well as abroad 
(TAES, US) that had been included in previous floristic inventories of the study area (Herrera 2001; Herrera & 
Peterson 2007; Herrera & Cortes 2009, 2010. Herrera et al. 2010; Cortes & Herrera 2011). 

f w 

Potential species distribution prediction was carried out using the Bioclim model included in the Di- 
va ~GlS program applied to Global Climate Data from the WorldCiim bioclimatic database (http://www. 
w orldclim.org/current). In brief, the Bioclim model identifies patterns of geographical distribution of species 
horn 19 bioclimatic variables (contained in the database) and determines areas where species could thrive. 


n like other models, this model uses presence-only data (Kenth & Carmel 2011, Rodriguez-Soto et al. 2011), 
characterizing the areas which are within the climate cover that have environmental characteristics similar to 
those occupied by the species under study. Bioclim can be used to describe the environment in which a certain 
species has been previously found in order to identify other areas where the species can live in the current 
condition as well as to identify areas where the same species could be present if the climate landscape alters. 
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llon of about one kilometer. It can be used to generate maps or create a spatial model in GIS or other computer 
software. This database was generated by the interpolation of observed and representative data from the years 

^0 to 2000 (Hij-mans et al. 2005). 


of prediction maps of di 


Potential distribution maps 



generated for the four grass 
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s Pedes with the Bioclim model, initially classifying the area into the following categories of probability of oc- 
cu rrence: Unsuitable, Low (0-15 percentiles). Medium (3.5-5 percentiles). High (5-10 percentiles). Very High 
( f0-20 percentiles), and Excellent (20-34 percentiles). In order to determine the main climatic characteristics 
rele vant to the species distribution, the maps were remodelled so that only cells/areas with high, very high and 

**®*®®ttprobability were shown. 













444 


Journal of the Botanical Research Institute of Texas 7(1) 


Creation of temperature and precipitation maps for each species .—Maps of the mean annual temperature 
and annual precipitation of the study area were generated from the climate database using Diva-GIS functions. 
Afterwards, using overlapping map functions, annual precipitation and mean annual temperature maps were 
generated for the high-probability distribution areas oi each species. 

Crossing of temperature and precipitation maps .—-The annual precipitation and mean annual temperature 
maps were cross-mapped so that for every species, each high-probability distribution map cell also showed the 
precipitation and temperature values. Furthermore, the cell frequency for the various precipitation and tem¬ 
perature values was obtained. 


RESULTS 


Maps of the predicted distribution of four grass species were generated with the Diva-GIS Bioclim algorithm. 
The selected species were Muhlenbergia peruviana, Muhlenbergia porteri, Bouteloua gracilis, and Melinis rcpcns. 
which were chosen because they have been collected in different ecosystems and our results show the varieu 
of climate conditions in which they can develop. Said maps show the level of probability of occurrence for each 
of the species in the study area (Guarino et al 2002; Hijmans et al. 2004). 

Muhlenbergia peruviana is listed as a species commonly distributed in pine-oak forest communities and 
the potential distribution map that was generated (Fig. 2) shows that areas of the Sierra Madre Occidental ha\I 
greater probability of occurrence for this species. 

When the mean annual temperature and total precipitation maps of high, very high and excellent prob¬ 
ability of distribution areas of M. peruviana were obtained, we found that this species is located in areas with a 
relatively low mean annual temperature between 12 and 17°C, particularly at 14°C. Furthermore, the annua 
precipitation ranged between 600 and 1000 mm, notwithstanding the fact that the largest area was comprise 
in places with 600 mm (Tables 1 and 2). With this in mind, it is possible to establish that the distribution of M 

peruviana is wit Irtn areas with low temperatures and relatively high precipitation. 

The potential distribution map of Bouteloua gracilis, which is known to be mainly distributed in natu 
grasslands and within some xerophilus scrubs (Fig. 3), shows a large area with high to excellent values of po¬ 


tential distribution. 


The distribution area of B. gracilis, with greater likelihood values, had annual precipitations ranging ^ 
400 to 800 mm, with an average of 560 mm although the most common precipitation values (mode) were J 
and 600 mm (these had the highest percentage of surface area). The mean annual temperature established in 
our model showed a range between 13 and 19°C, with an average of 16.2°C and a mode located at 17 C. 

The potential distribution map of Muhlenbergia porteri (Fig. 4), which is known to be present within .' 

rophilus scrubs (Fig. 4), placed this species within the driest regions of the study area. 

The area with the highest predicted distribution of M. porteri had a mean annual temperature rangt c 
between 12 and 21°C, with the mean located at 17.8°C and the mode at 18°C, while the annual precipkati 

ranged between 300 and 900 mm with a mean centered at 482 mm and a mode of400 mm. 

Climatic requirements for each of the three native species are summarized in Table 3, showing 
expected, they belong to different ecosystems. ^ 

The potential distribution map of the nonnative introduced grass Melinis repens (Fig. 5) showe 
had the smallest area of likely distribution. It is possible that these results are influenced by the fact tna 
species was represented in low numbers with all collections used in this study. 

Nevertheless, the distribution area with the highest likelihood of having Melinis repens showe ^ 

nual temperatures ranging between 14 and 18°C, with an average of 16°C and a mode also at 16°C. ^ 

nual precipitation reached an average of 609 mm, with a range between 500 to 800 mm and a mode o 


that, as 


DISCUSSION 


Results from previous inventories of grass species carried out in 

Mexico (Herrera 2001; Herrera & Peterson 2007; Cortes & Herrera 2010; Herrera & Cortes 
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Estate brat 


Colect of Muhlenberg* peruviana 
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Without data 




kilometers 


Fifi. 2. Map of the potential distribution of Muhlenbergia peruviana obtained with Biodim. 


rera et al. 2010; Herrera et al. in prep) allowed us to construct probable distribution maps of close to 500 species 
°f which the actual areas of distribution are known; this work shows the results from a subset of four selected 
species. The distribution of grass species in northern Mexico is mainly determined by climatic factors (Rze- 
dowski 1978); therefore it can be assumed that this same analysis applied to the remainder of the grass species 
within this area has a very high probability of increasing our knowledge on their potential distribution. 
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Table 1 . Proportion (%) of cells/areas with high probability of distribution of Muhlenbergia peruviana with regards to mean annual temperature and total precipitate 


Annual precipitation (mm) 
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12 

21 

28 

24 

13 

2 

0 

0 

0 

0 

100 


Table 2. Annual precipitation and mean temperatures associated with potential distribution of Muhlenbergia peruviana 



Num. of cells 


Minimum 


Maximum 


Mean 


Standardesv. 


Moda 



Annual precipitation 
Annual mean temperature 


2449 

12 


500 

17 


1000 

14 


733 

1.32 


142 

15 


600 



Table 3. Comparison of annual precipitation and mean temperature requirements of the three native species analyzed. 



Muhlenbergia peruviana 


Bouteloua gracilis 


Muhlenbergia porteri 



Annual mean precipitation 733 mm 560 mm 482 mm 

Annual mean temperature 14°C 16.2°C 17.8 a C 



In this study we modeled the potential distribution of four grass species from (he arid and semiand re 
gions in northern Mexico, allowing us to predict the area in which it is more likely to find each of them as 
as to determine the main climatic characteristics necessary for each species to be present. 

The maps of potential distribution of Muhlenbergia peruviana, Bouteloua gracilis, Muhlenbergia p°' u 
and Melinis repens, generated with the Bioclim model using 19 bioclimatic variables, show the probabi ity ^ 
distribution of these species within the study area in regard to their bioclimatic characteristics. These map 
lowed us to determine the annual precipitation and mean annual temperature conditions in which eachspe 
has been found. The possible distribution areas can be considered to be critical for further fieldwork as 1 

* A i ct 

show similarities in the main climatic characteristics with those areas where these species are a 

i ; ^aWnt cfr 

Our analysis clearly shows that the three native species that were selected grow under very u 
matic conditions, lending support to our initial decision of selecting species from different ecological**®^ 
Nevertheless, Melinis repens (the introduced species) was shown to share a large ecological niche wit 

1 r ^ * 

loua gracilis (the grasslands species) since the former has climatic requirements similar to those o 


such as 


forests 


This fact could mean that this African species would not be able to invade other communities ^ 

or xerophy tic scrubs, only posing a problem for grasslands. With our model it is possible to apprecia e 


tual distribution shared between a native species and an introduced species. 


th 


r _ . The hvp 0111 ' 

(The replacement of native species with less desirable introduced species affects biodiversity. 1 ; . 

esis that biodiversity has an influence on productivity and stability of the grassland ecosystem has 
ously tested in European grasslands by Hector et al. (1999) and in North American grasslands by 1 n ^ 
(1997). Both groups agree that the loss of biodiversity directly affects productivity and stability of g raSS 

i \ :q an inva* 

One of the grasses that was studied in our work, Melinis repens (Zizka rose Natal grass;, ntan in¬ 

species and the prediction of critical areas of future expansion is an important issue in order to P reve 
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Rc. i. Map of the potential distribution of Bouteloua gracilis obtained with Biodiin. 

vasion and later displacement of native grassland species in Mexico. Melinis repens is an African grass that was 
recently introduced into Mexico (less than 50 years ago), first found at the edge of roads, a fact that possibly 
helped this species to spread, and it now invades and displaces native species in grasslands in northern Mexico 
(Herrera Arrieta et al. 2011; Herrera Corral et al. 2011). This grass is highly aggressive and this characteristic 
could be allowing it to increase its cover and dominate the distribution areas of Bouteloua gracilis (blue grass). 
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Fig. 4. Map of the potential distribution of Muhlenbergia porteri obtained with Biodim. 

However, Figure 5 shows that climatic conditions could be containing Zizka rose Natal grass from sp 
into every space that blue grass now occupies within these grasslands (see Fig. 4). 


eadin- 


CONCLUSIONS 


1 1 * ' A- ^ 

With this study we have shown that it is possible to carry out prediction modeling of the potenti ^ 
or ecological niche of grass species based upon presence-only data (Kent & Carmel 2011; Rodriguez 
2011) that can be obtained from specimens collected with floristic inventory purposes. 
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he. S. Map of the potential distribution of Melinis repens obtained with Biodim. 


Also, this study demonstrates that the overlap in distribution between native and invasive spec.es can be 
shown, providing valuable information that can be used for preventing the spread of invasive species an e 
displacement of native species. We also uncovered evidence that it is possible that dhUOC requirements are 
preventing the expansion of the African-introduced species Melinis repens (Zixka rose Natal grass), although 

■as spread into grasslands in northern Mexico during the last 40 \ eat s. 
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Our work provides information on the potential distribution of grass species that could be useful in fu¬ 
ture biological assessments as well as for determining areas that are critical for future fieldwork. Furthermore, 
studies similar to this one should provide support for decisions relating to conservation and management of 
grasses in Mexico. 
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